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The Dirac fermion in the random chiral models is studied 
which includes the random gauge field model and the random 
hopping model. We focus on a connection between continuum 
and lattice models to give a clear perspective for the random 
chiral models. Two distinct structures of density of states 
(DoS) around zero energy, one is a power-law dependence on 
energy in the intermediate energy range and the other is a 
diverging one at zero energy, are revealed by an extensive nu- 
merical study for large systems up to 250 x 250. For the ran- 
dom hopping model, our finding of the diverging DoS within 
very narrow energy range reconciles previous inconsistencies 
between the lattice and the continuum models. 
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Dimensionalities and symmetries play central roles 
for universalities in the Anderson localization problem. 
More than two decades ago, Abrahams, Anderson, Lic- 
ciardello and Ramakrishnan |ll| presented the well-known 
scaling theory, predicting that electron wave functions 
always localize in one and two dimensions and metal- 
insulator transition occurs in three dimensions. How- 
ever, even in two dimensions, delocalized states are also 
marginally allowed to appear when systems possesses 
some special symmetries. For symmetries, it is often con- 
venient to borrow terminologies from the random matrix 
theory which was first introduced by Wigner and Dyson. 
Recently, Altland and Zirnbauer g] have reported seven 
new symmetry classes in conjunction with mathemati- 
cal classification scheme of the Riemannian symmetric 
spaces. 

Amongst these new symmetry classes, chiral models 
have attracted much attention as a novel exception for 
the scaling theory. Compared to the conventional mod- 
els where randomness enters as an on-site potential, it 
resides on links (i.e., as a gauge field) in chiral models. 
This type of randomness may play important roles for the 
composite fermion theory of fractional quantum Hall ef- 
fects and vortex states of dirty superconductors. In view 
of a localization problem, models with this randomness 
have a special symmetry, referred as chiral symmetry, 
and thereby belong to a new universality class B. This 
symmetry is expected to affects localization properties of 
the systems drastically. 

Several chiral models on a lattice, which are convenient 
for numerics, have been studied and interesting physics 
around zero-energy have been revealed. DoS of these 
models show singularities at zero-energy and the corre- 
sponding wave functions exhibit a delocalized multifrac- 
tal behavior. Examples of these include the Gade's model 



, the random flux model [^ , and the tt-Aux model with 
link disorders [p|. These models, defined on a two dimen- 
sional square lattice, have the chiral symmetry which is 
conveniently stated as {H, 7} = where 7 is a matrix 
which changes the sign of wave functions on a sub lat- 
tice. Consequently, for any realization of the disorder, 
the energy spectrum is invariant under the transforma- 
tion E —^ —E. Therefore, given an eigenstate ip with 
energy E, jip is also an eigenstate with energy —E. This 
symmetry is responsible for the existence of delocalized 
states at zero-energy. 

Chiral models in a continuum space have also been 
investigated extensively. Especially, models which in- 
clude the Dirac fermions have attracted an interest. The 
Dirac fermion is a quasi particle which appears in sev- 
eral interesting situations in condensed matter physics 
such as d-wave superconductivity, graphite sheets, the 
gap-closing transition in quantum Hall effects Q, the 
Chalker-Coddington network model J7|, and the mean 
field theory of the t — J model pj. Effects of random- 
ness are also of fundamental interest in these contexts. 
Since many analytic approaches such as a field theoretic 
one are applicable for the Dirac fermion, several inter- 
esting results have been obtained so far. Amongst them, 
there exist simple models which allow us to construct an 
explicit zero-mode wave function for any realization of 
disorder ^. Due to this advantage, it has been revealed 
that the zero-mode wave functions are not localized and 
exhibit a multifractal behavior like ones in the lattice 
models. Moreover, there exists a transition in the multi- 
fractal spectrum pX)| and the density-density correlation 
||ll| as the disorder strength varied. 

Although the delocalized multifractal nature of the ex- 
act zero-energy states has been well established now, our 
knowledge for finite energy states, especially for DoS 
around zero-energy, is still in confusion. For example, 
a continuum model where the Dirac fermions feel ran- 
dom gauge fields was studied in Ref. 1^,0 and it was 
found that DoS exhibits a power-law dependence on en- 
ergy. On the other hand, a similar model with species 
doubling was studied in Ref. [Of and diverging DoS was 
found. From the lattice point of view, the corresponding 
model including Dirac fermions has not been studied so 
far. However, the random flux model, a lattice model 
where flux is randomly distributed for each plaquette, 
shows diverging DoS ||l^. A clear link between them is 
missing. 

Another example is the case where the Dirac fermion 
feels imaginary vector potentials. For this case, there 
even exist inconsistent results between lattice and contin- 



uum models. Field theoretic studies ||1^,|15[ predicts di- 
verging DoS at zero-energy for any randomness strength. 
However, a numerical result for the random tt-Aux model, 
where a random hopping amplitude act as a imaginary 
vector potential, shows that DoS behaves as a power-law 
with its exponent dependent on the disorder strength for 
weak randomness H . 

In this paper, we make an attempt to clarify the re- 
lationship between lattice models and continuum Dirac 
fermions, reconcile some inconsistent results in the pre- 
vious studies, and thereby give a whole perspective for 
the localization problem of the chiral models. For these 
purposes, we consider lattice models which recover the 
random Dirac fermions in the continuum limit. We focus 
on DoS, especially around zero energy where quantum in- 
terference play an important role and fully quantum me- 
chanical treatments are necessary. We use the transfer- 
matrix method developed in Ref. ||l^. It allows us to 
handle large enough systems up to 250 x 250 , which 
is indispensable to give reliable results for localization 
problems. 

We realize Dirac fermions on a 2D square lattice via 
the TT-flux model: 
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(-)^" , t'j+'-g,] = 1 ^'^'i ^^^ piercing 
a plaquette is it. The energy spectrum is given by 
E = ±2-^cos2 k^a 



cos^ kyQ where k belongs to the 



magnetic Brillouin zone (— 7r/a,7r/a] x (0,7r/a] and a 
is a lattice constant. In the continuum limit a -^ 0, 
this includes the doubled massless Dirac fermions around 
k = (7r/2a, 7r/2a) and (— 7r/2a,7r/2a). This realization is 
a minimum model for our purposes. 

In the following, we will consider two types of disorder 
which live on a link: the random gauge field and the 
random hopping. We implement the random gauge field 
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Taking a Coulomb gauge, we determine the random 
gauge field Aij via a scalar potential $ on a dual lat- 
tice as Aj+i^j = ($j+(£+j;)/2 - ^] + (x-y)/2) /a.Aj+y^j = 

{<^j+{x+y)/2 - ^j-{x-y)/2) /a where $ is randomly 



chosen from a simple Gaussian distribution P['^] oc 
exp[— |- X)<i7> i^i ^ ^j) l'^\- III the continmmr limit, 
this reduces to P[<I>] oc exp[— ^/d^a;(V$) ], a natural 
choice for effective field theoretic treatments. As stated 
above, this Hamiltonian possesses the chiral symmetry 
{H,7} ~ due to the special nature of randomness. 
Note also that time-reversal invariance is broken by the 
random magnetic flux. 

In the continuum limit, the Hanriltonian is expressed 
as 1,0 
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where D := —idxCTx — idyay + Ax{yi)ax + Ay{x.)ay + 
M{x.)az + ^(x)l , cTi=x,y,z are the 2x2 Pauli matri- 
ces, and ^ is a four-component spinor. The coefficients 
Ax, Ay, M and V are arbitrary complex fields. This is the 
most general chiral symmetric form of a Dirac Hamilto- 
nian in a continuum space. The chiral symmetry is, in 
the present basis, expressed as {Ti, cr^ (8) 1} = 0. For the 
random gauge field model, {Ax, Ay) is real and serves as 
the random gauge field. The other coefficients M and V 
do not appear at the first order in the lattice constant a, 
but they are non-zero for the lattice models in general. 
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FIG. 1. Density of states for the vr-flux model with a 
random gauge field on a 50 x 50 lattice for g — 0.2 — 1.6 (from 
bottom to top at _E = 0). The small imaginary part of energy 
S is 0.01. Quenched averaging is taken over 30 samples. 
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FIG. 2. Same as Fig. |l[ on a 100 x 100 lattice for 
g = 0.1 — 1.8, and 5 = 0.005, averaged over 40 samples. Inset; 
same as Fig. on a 250 x 250 lattice for g = 0.4, 0.45, 0.5, 
and 5 — 0.0005, averaged over 50 samples. 

For strong enough randomness, this model is naturally 
related to another lattice model, the random flux model. 
(By the random flux model, we mean the model discussed 
in Ref. 0,0 which does not include Dirac fermions, while 
we simply call the model with Dirac fermions defined 
above as the random gauge field model.) To see this, 
note that the spatial correlation of the flux per plaquette 
4>i := a^Ai$i is expressed as {4'i4>j) = a^AiAjG{\i — j\) 
where G{\i — j\) := {^i^j) is a lattice Green's function 
(2) for <&, and A^ is lattice Laplacian at a site i. Taking 



a continuum limit a — > 0, we find {4>i4>j) — for finite 
a|i — j'l, i.e., the spatial correlation of flux is short ranged. 
Moreover, for strong enough randomness, lattice fermions 
do not remember original flux tt anymore, and we expect 
this model to exhibit a similar behavior to that of the 
random flux model where the flux for each plaquette is 
an independent random variable. 

An effective sigma model via supersymmetry (SUSY) 
technique was studied in Ref. |1^, predicting di- 
verging DoS for all range of randomness strength. 
A similar model without species doubling Ti = 
J d^a;^^ [cr • p + cr • A] ^ has also been studied [^ 12 1 and 
power-law DoS with a disorder-dependent exponent was 
proposed. On the other hand, the random flux model, 
which is expected to be the strong randomness limit of 
the present case, shows a diverging DoS MS. 
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FIG. 3. Density of states for the random hopping model 
on a 50 X 50 lattice for g = 0.3 — 0.7 (from bottom to top 
at _E = 0). The small imaginary part of energy S is 0.01. 
Quenched averaging is taken over 30 samples. 
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FIG. 4. Same as Fig. g on a 100 x 100 lattice for 
g = 0.3 — 0.7, and 5 = 0.005, averaged over 50 samples. Inset: 
same as Fig. | on a 250 x 250 lattice for g = 0.5,0.51,0.52, 
and 5 = 0.0005, averaged over 70 samples. 

We also consider the case where fermions on a 2D 
square lattice with 7r-flux per plaquette feel random hop- 
ping amplitudes: 



,rh _ .pure 



St, 



(3) 



where Stij is a real random variable. As for a probability 
distribution of St, we consider a Gaussian distribution 
P[St](xexp[-{5t)y2g^]. 

In this case, {Ax,Ay) in Eq.(y) is purely imaginary, 
and Dirac fermions feel an imaginary random gauge po- 
tential. Note that the total four-component Hamiltonian 
is hermitian. Unlike the random gauge field model, time- 
reversal symmetry is not broken for this model since we 
can take all matrix elements to be real. 

As is the case of the random gauge field model, this 
model is connected to another lattice model for suffi- 
ciently strong randomness. For strong enough random- 
ness, again, one expects that fermions do not remember 
the original flux anymore and shows a similar behavior 
to that of the Gade's model. 

Fukui [n5| studied the continuum model by a replica 
nonlinear sigma model with a large number of fermion fla- 
vors and proposed that DoS at zero-energy diverges and 
delocalized states exist at the band center for any ran- 
domness strength. One-loop renormalization group (RG) 
study with SUSY method was also applied for this model 
in Ref. |13], predicting diverging DoS at zero energy as 
p{E) - ^^-c^HEr/e) ^Yiere Er is a constant and c 
depends on randomness. A conjectural RG flow beyond 
one-loop order was presented in Ref. |17]| and similar re- 
sults were obtained. The Gade's model, which we expect 
to be the strong randomness limit of the present model, 
also shows diverging DoS [|j. On the other hand, accord- 
ing to a numerical study for the corresponding lattice 
model, vanishing DoS p{E) ~ \E\°' was observed with a 
disorder-dependent exponent a > for weak randomness 
H . We re-examined this problem by the transfer-matrix 
method. 

A quantity of interest we investigate numerically is ran- 
dom averaged DoS: {p{E)) := {J;IY.^S{E^ - E)). For 
the pure tt-Sux model, DoS vanishes linearly around 
zero-energy due to the relativistic dispersion of Dirac 
fermions. Since for the strong enough disorder, the 
present models are expected to exhibit similar proper- 
ties to those of the random flux model or the Gade's 
model, both of which exhibit divergent DoS at zero- 
energy 1^,0, a natural question is how the vanishing 
DoS at zero-energy becomes divergent as we increase the 
disorder strength g. It should be contrasted to the usual 
case where DoS of pure systems are also divergent at 
zero-energy due to the van-Hove singularities of the tight 
binding model on a two-dimensional square lattice. 

First, we discuss the random gauge field model, i.e., 
TT-flux model with a random gauge field as in Eq. (|lj). 
In Fig. ^ and 0, we present numerically calculated DoS 
for several g. In the transfer matrix method, the Green's 
function for a given energy E is calculated at £^ -I- iS. We 
chose S ^ 0.01 — 0.0005 which gives us enough resolution 
of energy. We used a square geometry rather than a 
quasi one-dimensional one since we are mainly interested 
in pure two-dimensional properties. 

As shown, there exist two structures. For over the 



wide range of the energy scale, the global feature of DoS 
changes from V-shaped to flat one, the latter is character- 
istic to the random flux model as expected. In addition, 
if we look more precisely, for very tiny region around 
zero-energy, there exists another structure, i.e., diverg- 
ing DoS for sufficiently large g. Note that since we used 
a lattice with even number of sites and adopted periodic 
boundary condition, there is no exact zero energy state, 
i.e., the diverging DoS found here is not an artifact by 
the special choice of boundary conditions. The numerical 
results found here are consistent with the sigma model 
study via SUSY technique in Ref. [|3|, where the diver- 
gence of DoS is predicted. 

Now let us go on to the tt-Aux model with random 
hopping amplitudes ( Eq.(|^) ). Calculated DoS for sys- 
tems up to 250 X 250 are shown in Fig.0 and 0. Again, 
one recognizes two different structures as in the random 
gauge field model. Away from zero energy, DoS exhibits 
a power law behavior with a disorder-dependent expo- 
nent as in Ref. H . The whole DoS profile becomes diver- 
gent as the disorder strength increased, which is natural 
since the present model is expected to show similar be- 
haviors to that of the Gade's model. On the other hand, 
as is the case with the random gauge field model stud- 
ied above, within a very narrow range near zero-energy, 
there exists another structure. As the disorder strength 
increased, DoS divergence at zero-energy occurs when g 
reaches about the order of the bandwidth. This behavior 
is independent of the global power-law profile away from 
the zero energy and before it turns to be divergent. These 
findings are consistent with the field theoretic analyses, 
especially with the SUSY approach where the divergence 
of DoS is predicted and a power-law behavior appears in 
the intermediate region of the RG fiow. We also investi- 
gated the case where 5t is uniformly distributed between 
[— w/2,w/2]. The results for this case are qualitatively 
similar to those of a Gaussian distribution. 

We should compare the present case to dirty d-wave 
superconductors which also include Dirac fermions, how- 
ever, the underlying symmetry class is different in the 
random matrix theory yj. For dirty d-wave supercon- 
ductors, all quasiparticle states are found to be localized, 
and a small energy scale around zero energy appearers 
where quantum interference effects produce several criti- 
cal DoS profiles depending on the details of randomness 
p8[ . This energy scale is determined by the diffusion con- 
stant in the diffusive regime and the localization length 
of quasiparticles [|l^. For the present case, on the other 
hand, zero-energy states are delocalized due to the chi- 
ral symmetry, and quantum interference effects give rise 
to diverging DoS at zero energy. The emergence of a 
power-law behavior away from zero energy is, from the 
effective field theory point of view, well described by the 
intermediate regime of the RG fiow |1^ . 

In conclusion, we have investigated DoS of the chiral 
models where the chiral symmetry plays crucial roles for 
the existent of the delocalized states at zero energy. We 
considered the lattice counter parts of continuum mod- 



els including the Dirac fermions, the random gauge field 
model and the random hopping model, which have been 
well studied by field theoretic methods. Large-scale cal- 
culations by the transfer matrix method up to 250 x 250 
lattices revealed the existence of two distinct structures 
in DoS, one is a power-law behavior in the intermediate 
energy range and the other is a diverging one around zero 
energy. Quantum interference plays crucial roles for these 
fine structures, which naive semiclassical treatments may 
miss. Our finding of the diverging DoS at zero-energy for 
the random hopping model reconciles inconsistencies be- 
tween lattice and continuum models and thus establishes 
a clear connection between them. 

We thank Y. Morita for fruitful discussions. The com- 
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Computing Facility and at the Supercomputing Center, 
ISSP, University of Tokyo. 
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